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By investigation of the dielectric response of a Fermi-Dirac plasma in the linear limit and eval- 
uation of the electrostatic potential around the positive stationary test charge, we find that the 
Shukla-Eliasson attractive force is present for the plasma density range expected in the interi- 
ors of large planets for a wide range of plasma atomic-number. This research which is based on 
the generalized electron Fermi-momentum further confirms the existence of the newly discovered 
Lennard- Jones-like attractive potential and its inevitable role in plasma crystallization in the cores 
of planets. Moreover, it is observed that the characteristics of the attractive potential is strongly 
sensitive to the variation of the plasma density and composition. Current research can also have 
applications in the study of strong laser-matter interactions and inertially confined plasmas. 
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INTRODUCTION 



Miniaturization in technology on one hand and interest in astrophysical phenomena on the other have motivated a 
great deal of enthusiasm towards the study of quantum plasmas. In a series of papers Fowler, Chandrasekhar, Bohm, 
Pines, and Levine (@, EI H H3) have founded the concepts of relativistic-degeneracy and quantum plasma 
theories. Recent developments of quantum hydrodynamics (QHD) ([Ullifij]) and spin-1/2 quantum magnetohydro- 
dynamics (QMHD) models ( [1, [H HH ) based on Wigner-Poisson and Schrodinger-Poisson formulation (e.g. see[l9j]) 
have opened a new area of research in the field of cold ionized plasmas. Quantum hydrodynamic theories can be 
used to evaluate the dynamics of spin-1/2 electron- hole systems such as semiconductors ([H,[2^]) and nanostructured 
materials ([Hj]). Recent studies based on the QHD model have revealed that the particle nonlocality feature intro- 
duced into the quantum plasma model due to the Bohm force is a new feature only present in Fermi fluid models 
and can introduce additional critical quantum effects in the collective dynamics of Fermi-Dirac plasmas. Quantum 
plasmas have much reduced average interparticle distances compared to that in the classical counterparts. It has 
also been shown that the particle diffraction effect in the quantum plasma can lead to new multistream plasma in- 
stabilities not encountered in ordinary classical plasmas (plj)- Extended quantum hydrodynamics formalisms have 
recently appeared for magnetized and unmagnetized cases taking into account the relativistic spin-orbit interaction 
of Fermi-Dirac particles ubiquitous in strong laser-matter interactions and astrophysical dense plasmas ([H l20l l32l|). 



Applications of QMHD in quantum plasmas, on the other hands, reveals that dynamics of the Fermi spin-1/2 
plasma is fundamentally affected by the pressure caused due to the electron spin and orbit magnetization ( 10, 30|V 
Furthermore, it has been shown that the magnetization in quantum plasmas can give rise to distinct degeneracy regimes 
(0) and remarkable behavior for spin-induced nonlinear waves (Q). One of the distinguished manifestation of electron 
quantum diffraction effect in a degenerate plasmas is the plasma dielectric response to the external electromagnetic 
perturbations. The plasma dielectric susceptibility for unmagnetized and magnetized quantum plasmas, have recently 
been investigated by many authors. The new findings indicate the critical influence of quantum effects on the Debye 
shielding in Fermi plasmas [36[. as a fundamental problem, there has been many previous investigations on the 
potential around a test charge in variety of plasma environments ( (3714391 EH). More recently, by evaluating the 
electrodynamic potential of a moving test charge in a quantum plasma a modified Debye-Hiickel (DH) potential, i.e., 
the ordinary DH potential coupled with an oscillatory term have found in quantum plasmas (40| . Such modified DH 
potential, particularly caused by the presence of electron degeneracy, interaction energies and the electron nonlocality 
effect, has been shown to lead to a new type of Lennard- Jones-like attractive force with a pronounced potential valley 
near the test particle ((4ll|). The discovery of such new force in quantum plasmas has fundamental significance for the 
plasma coupling, crystallization and quantum transport effects in dense plasmas. The introduction of such attractive 
force, which brings the ions closer together, may even lead to breakthroughs in semiconductor, supercomputing and 
nanotechnology sciences. 



The astrophysical applications of the quantum hydrostatic model became apparent, since the pioneering work of 
Chandrasekhar ( |lll4l3j ) on the relativistic degeneracy. Due to the change in the equation of state of the electron 
degeneracy pressure in the relativistic degeneracy limit a fundamental limit on white dwarf-mass has been discovered 
( 14|) and distinct features has been shown to exist in the QHD wave dynamics of ultrarelativistic Fermi-Dirac plasmas 
( 4j). On the other hand, the application of QMHD model suggests some large-scale effects such as modification on 
the well-known Jeans criteria [25| due to the presence of electron spin and orbital magnetizations. Therefore, it is 
appealing to examine the effect of relativistic degeneracy on quantum plasma shielding and its dielectric response. 
In a relativistically degenerate Fermi-Dirac plasma Salpeter ([42j) has evaluated the effects of variety of interaction 
potentials such as Coulomb lattice-energy, Thomas-Fermi distribution, electron-exchange, and correlation on stellar 
hydrostatic stability for a wide range of relativity parameter and has found more significant contribution due to 
Coulomb lattice energy compared to others interaction effects. Such a plasma interaction effect although being a 
few percent of the dominant electron degeneracy effect, can give rise to crystallization of the Fermi-Dirac plasma in 
superdense stellar cores. In this paper we aim at examining the linear plasma response and the potential of a test 
charge due to the presence of Coulomb interactions and the Bohm force in relativistically degenerate Fermi-Dirac 
plasmas for a wide range of relativity parameter (electron number-density). We present the quantum hydrodynamics 
(QHD) model in Sec. „ and examine the plasma dielectric constant in the linear limit in Sec. „ Then, we calculate the 
modified Debye-Huckel potential around a stationary positive test charge in Sec. „ giving the parametric description 
of the associated potential in Sec. D Finally, we make numerical evaluation of the calculated potential based on 
independent plasma parameters in Sec. n and draw our conclusions in Sec. „ 
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QUANTUM HYDRODYNAMICS MODEL 



In this model we consider a plasma of nondegenerate stationary ions with charge Z quasineutrally distributed into 
a completely degenerate Fermi-Dirac electron fluid. The ion thermal energy is neglected compared to the electron 
Fermi-energy (Ep e ), hence, they are assumed to be static and unperturbed to external simulations compared to light 
electrons (m e /mi -C I). Furthermore, we use the nearly-free electron model neglecting the electron fluid viscosity. 
The electron viscosity has been shown to becomes dominant in much larger densities such as in neutron-star crusts 



2l| . On the other hand, it is also that in quantum plasmas the electron-electron collisions are limited due to the 
Pauli blocking-mechanism ([35]). In the electron number-density regime, considered in this work (I0 22 — 10 26 /cm 3 ), 
one may also assume infinite plasma conductivity (a ~ 10 22 ) thereby neglecting the ion-electron collision frequency 
([33:j) compared with the characteristic plasma frequencies. Therefore, the hydrodynamic evolution of Fermi-Dirac 
electron gas is governed by the following set of equations (0), as follows 

% + V ■ (mi) = 0, 

f + (u • V)u = eV0 - ^VP G + ^ V [A^i/V^ , (1) 
Acj)(r) = 4Tre(n — n ), 

where, fi and c are scaled Plank-constant and speed of light in vacuum, respectively. Also, Pq is the generalized 
pressure of plasma composed of electron degeneracy, exchange and Coulomb interaction pressures, i.e., Pq = Pd + 
Pec + Pc an d no is the unperturbed electron density. Ignoring the negligible ion pressure compared to that of the 
quantum pressure of electrons, a general form of electron degeneracy pressure for a wide range of plasma density is 
given as ([12j) 



P,. = 



3/i 3 



C — \R (2R 2 - 3) y/l + B? + 3 In [p + yj\ + i? 2 ] } , (2) 



where, the Chandrasekhar relativity parameter, P, is related to the Fermi relativistic-momentum via the relation; 
R = PFe/rrieC — (n/ric) 1 / 3 with n c ~ 5.9 x 10 29 cto~ 3 being the normalizing value. Furthermore, a relativistic 
degeneracy parameter, Rq, may be introduced as R — Ro(n/no) 1 ' 3 where Ro — (no/nj 1 / 3 with uq being the 
unperturbed equilibrium electron number-density. It can be easily shown that, for the degenerate electron pressure 
represented as a polytropic form of P = P^n/no) 1 the values of 7 = {5/3,4/3} correspond to the degeneracy limits 
R = {0, 00} of Eq. ([2]) as nonrelativistic and ultrarelativistic degeneracy limits, respectively. The change in equation 
of state due to the relativistic degeneracy of electrons has been shown to lead to the well-known white-dwarf mass- 
limit (13]) and distinct features in nonlinear wave dynamics in Fermi-Dirac plasmas (0). It can be observed that 
the polytropic form of the Chanrasekhar pressure approaches nicely to the standard Fermi degeneracy pressure for 
the extremely nonrelativistic number-density range [5j. It is evident that one needs to use the full form of pressure, 
valid for a wide plasma density-range, in order to evaluate the Shukla-Eliasson force and phase separation in dense 
quantum plasma to be shown present also in the interiors of large planets. 

As the plasma becomes denser, the corrections to the equation of state of the plasma due to the particle interactions 
become important. Such interaction terms in the strong coupling plasma limit are due to the electrostatic plasma 
interactions (Coulomb interactions), electron exchange, electron nonuniform distribution (Thomas- Fermi correction), 
ion correlation and some other minor effects ((42|). It is known that the plasma coupling parameter, T, is inversely 
proportional to the interparticle distances, r^, and the temperature T of the plasma species. Hence, with the inter- 
particle distances being of the same order for ions and electrons in the quasineutrality limit, the coupling parameter 
for ions is observed to be much larger than that for degenerate electrons, since, the interacting electron temperature is 
of the order of the Fermi-temperature being much higher than the ion thermal temperature (I\ ^> T e ). The Coulomb 
attractive potential, being the dominant interaction potential, is the main cause of the crystallization of Fermi-Dirac 
plasmas, under strong coupling conditions. Its dominant effects is well-understood in ordinary solids-state matter 
and is believed to play a fundamental role in Fermi-Dirac plasmas lattice formation. Calculation of the collective 
electrostatic energy for crystals with different periodic configurations such as the fee and the hep lattice structures 



confirm that the Coulomb energy of a lattice is minimally related to the lattice structure ([22j]). Salpeter (|42j), using 
a spherical noninteracting Wigner-Seitz cell of radius, r (l/n — 47T7o/3), has given an expression for the Coulomb 
energy of a relativistically degenerate plasma with ions forming a quasi-lattice at very high electron-density plasma 
regime such as in a white dwarf core or neutron star crust. It has been shown that such attractive force persists and 
even becomes stronger as the value of the relativity parameter or the electron number-density increases. On the other 
hand, the plasma coupling parameter which is given by T = (Ze) 2 / (roksT) may also be approximated in terms of 
the relativity parameter as; T ~ 0.23Z 5 / 3 P/Tg, (where T$ is the plasma temperature in units of 10 8 degrees Kelvin) 
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([2 If). It is known that, the plasma is in liquid form for the coupling parameter range V < 178 and in solid phase 
beyond this value ( 21]). Hence, it is concluded that for a typical planet core temperature, say 24000A" for Jupiter 
core, and a typical relativity parameter of R ~ 0.1, the coupling parameter value exceeds the liquid critical limiting 
value. Therefore, it is observed that, the planet cores for a wide range of plasma atomic-number value would be in 
a quasi-lattice phase. The effects of interaction potential on ion wave dynamics in Fermi-Dirac plasmas has been 
investigated recently @|. 



Pc 



ft 3 



IOtt 2 



1/3 



R 4 



(3) 



where, Z is the atomic number of ions and h is the Planck constant. It should be noted that, the Coulomb pressure 
presented here only gives a qualitatively right order for Z 2//3 (ro/rs) ~ 1 in the nonrelativistic plasma density limit 
(10 25 < uq/ cm 3 < 10 27 ), where, tq and tb are the Wigner-Seitz and Bohr radiuses, respectively. Therefore, in the 



proceeding we will confine our analysis to this validity limit which is relevant for the planet core density range, 
the other hand, the contribution from the electron exchange may be written as, 



P 



{h (Z? 4 + + I 32 + - I (P 2 - r 2 ) ln/3 + § (In/?) 2 



On 



(4) 



where, f3 = R + y/l + R 2 . The Coulomb interaction pressure is always negative, while, the pressure due to electron 
exchange changes the sign at the relativity parameter value around R ~ 2 [6(. It should be noted that in the 
forthcoming calculations we have included the generalized form of the exchange energy, given by Salpetcr, instead of 
that in (|4l|). since, the generalized exchange energy is valid for a wide electron number-density range. We further 
have neglected the effect of electron-correlation effect due to the vanishingly small contribution in the considered 
density regime (|42T|). 



LINEAR DIELECTRIC RESPONSE 



In order to obtain the dielectric constant for the present plasma model we use a standard procedure employed by 
Shukla & Eliasson ([41]). The QHD equations given in Eqs. (j6|) now read as 



at 
at 



V • nu = 0, 

(u • V)u = ^-V</> - c 2 V [* d + * xc + * c ] 
A(f>(r) — Ane(n — no), 



h 2 y A^/n 
2m'j, -Jn ' 



where, "f's arc the effective plasma potentials defined via the relation 



rrif, 



I? 3 



m, 



R 3 dR 



The plasma effective potentials are, then, given as 

«o 



l + i? 2 n 2 / 3 , t> xc = 



2tt 



Ron 



1/3 



3sinh _1 (i?o« 1/3 ) 
y/1 + R^n 2 ^~ 



-j3Ron 



1/3 



with 



_ 2a / 2Z \ 

P 5 



2/3 



do 



1/137. 



(5) 



(0) 



(7) 



(8) 



To evaluate the linear response of the plasma to harmonic excitations of form cx exp[i(k • r — wt)], we use the well- 
known harmonic operators V = ik and dt = —ioj to linearize the equation set Eqs. Assuming a small density 
perturbation around the equilibrium state (i.e. u = ui, <fi = (pi, and n = no + n\ with n\ <C no), we obtain the 
following equations 



k • ui = 0, 



umi 

WUi — 



h k 2 kn x - c 2 Tk ni 



T = 



"0 

2tt 



3 



flgsinh" 1 .^ 
(l + Hg) 3 ' 2 



(9) 
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the reduction of which yields the first order potential perturbation and the electron dielectric susceptibility, as 



Xe 



Ami 



(10) 



where, \e is the dielectric susceptibility related to the dielectric constant as e(u), k) = 1 + Xe- I n the quasistationary 
limit, we obtain the normalized dielectric constant, as ((4l|) 



e(0,k) = 1 + w 



Ami 



k A + c 2 T 



(11) 



where, oj pe = \J AirriQe 2 /m e is the electron plasmon frequency. Note that, in many recent quantum plasma literature 
the use of the quantum diffraction parameter, H — h 2 uj 2 i /(2m e c 2 ) called the scaled ion-plasmon frequency, has become 

3/2 

widespread. The value of this parameter can be approximated as H ~ 0.0044Zi? which for metallic iron density of 
~ 10 23 /cm 3 and higher density of our concern (10 26 /cra 3 ) with same composition take the values of H ~ 0.00005 and 
H ~ 0.002, respectively. However, the corresponding values are much smaller than that suggested previously (18]) 
for metallic solids which is due to the normalization factor used in this work. 



PHASE SEPARATION IN FERMI-DIRAC PLASMAS 



Let us consider a positive stationary test charge in the plasma with Fermi-Dirac degenerate electrons. The shielding 
of such a test charge is a fundamental property of the plasma and can give useful information about the bound states 
between ions. The Debye-shiclding or the near-field electrostatic potential usually appears as 4>(r) oc (l/r)exp(— r/rd) 
where is the Debye-length and r is the distance from the test charge ([3]). In the extreme nonrelativistic degeneracy 
limit the plasma shielding potential can attribute to also a far-field inverse cubic dependence (r -3 ) from an outside 
inertial observer point of view (0)- L et us f° r simplicity assume a stationary test charge Qt = Q located at the 
origin. It is obvious that analysis can be extended to a moving test charge with velocity vt relative to the laboratory 
frame (|4p|V The Poisson equation in the presence of point-like stationary positive test charge can be written as, 



A<j>(r) = Ane(n - n Q ) - A-kQ6(y), 



(12) 



where, 8(r) is the Dirac delta- function ensuring the localization of the test charge at the origin, i.e., 8(0) = 1 and 
8(r ^ 0) = 0. Using the Fourier analysis in space, the electrostatic potential around the test charge can be rewritten 
in the form 



Q f exp(ik • r) 
k 2 e(k) 



2tt 2 



d A k, 



(13) 



where, e(k) is the plasma dielectric constant, given by Eq. (jlOp . We then appropriately decompose the integral 
following the procedure given by [4l[. Therefore, the integral in Eq. (fT3"|) can be decomposed into different parts, as 
follows 



[ ' Air 2 



where, the parameters fci.2 and L are given as 



l + L 



1 -L 



k ~\~ Aj-^ k ~\~ k^ 



exp(ik • r)d 3 k, 



,2 _ /2 lTv / l-4a 



L = Vl - Aa. 



with the parameters ko and a given as 



h 2 u;i 



k = 



Am 2 e c 4 T 2 ' '" u cVT 
Now, using the condition, <f>(r — oo) = 0, and the standard definition 



exp(ik • r) 



d A k = 



2ir 2 



(14) 



(15) 



(16) 



(17) 
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we may write the general solution for the first-order potential perturbation around the test charge as, 

4>{r) = f [(1 + L) e-^ r + (1 - L) e~ k ^] . (18) 

We, then, are in a position to evaluate the potential around the test charge in the Fermi-Dirac degenerate plasma 
with respect to two independent plasma parameters, namely, the relativistic degeneracy parameter, i?0i representing 
the normalized number-density of electrons and the atomic-number, Z, of the plasma ions. As it is evident, this can 
be only done in a numerical scheme presented below. We will present the phase separation in Fermi-Dirac plasma 
for different values of the plasma atomic-number using the scaled Bohr-radius, r s = Tb/tq, where, rs — Ti 2 /m e e 2 
is the Bohr-radius. This implies the relation; i?o ~ 0.014r s between the degeneracy parameter, Ro, and the scaled 
Bohr-radius, r s . 



SHUKLA-ELIASSON ATTRACTIVE POTENTIAL IN FERMI-DIRAC PLASMAS 



It is evident that for a physical solution to exist, in Eq. ([TBI , the value of the potential must be real. It can be easily 
confirmed that for all values of a > 1/4 the parameter L is imaginary. Therefore, for a > 1/4 (it can be confirmed that 
the parameter a is always positive), we have = 5ft(fc2) and Sy(fci) = — ^(fe), hence, — K(fci) ± 3(fci) ( t 41]) 

and the electrostatic potential can be written in the form <fi(r) = (Q/r) exp[— 5R(fc2)r]{cos[5R(/c2)?'] + $t(L) sin [^(fe )?*]}, 
where -ft and 3 refer to the real and imaginary parts of the parameter ki- This attractive Lennard- Jones-like potential 
relates to the ion bound-states of the Fermi-Dirac plasma. However, for a < 1/4, where the parameter L becomes 
real-valued, the potential of a test charge reduces to the repulsive Debye-Huckel or Yukawa-type potential. It is 
important to note that the Shukla-Eliasson force is an elegant manifestation of pure quantum mechanical effect and 
a consequence of the electron collective quantum behavior, since, it disappears in the h — > limit. It is also observed 
that, the characteristics of such attractive force depends strongly on the atomic-number of the plasma ions. As it will 
be revealed in Sec. „ our research highlights the importance of quantum tunneling effect on electrons through the high 
plasma number-density and the Coulomb negative pressure. It is also observed from Eq. (|16p that below a critical 
electron number-density, n$ cri for a given atomic-number, Z, the parameter T can becomes negative leading to an 
imaginary value for k$ and consequently for the screening potential. Therefore, there is a lower plasma number-density 
bound, where the attractive force can exist. This is completely natural, since, as the plasma density is lowered the 
degeneracy limit of Fermi-Dirac plasma is reached beyond which the attractive force should disappear. 

It is most surprising, however, to see some recent literature ( 45]) to question such rigorously furnished theory as 
quantum plasma theory and its fundamental applications, based on fallacious arguments. Also, based on the density 
functional theory computer simulations some authors (|8j) have argued against the applicability of linearized quantum 
hydrodynamic theory and the existence of Shukla-Eliasson attractive potential. One must note that, the development 
of improved functional theory is currently a very active area of research and it is far from clear whether the research 
will be successful in providing the substantial increase in accuracy desired ( (43| ) . In the following numerical analysis it 
is confirmed that such attractive potential not only persists in high density matter it also may be one of the essential 
ingredients causing the strong ion-coupling in the cores of planets. 



NUMERICAL ANALYSIS OF PHASE-SEPARATION 



Figure 1 shows the phase separation parameter, a, for the extended quantum plasma model in terms of the scaled 
fermion spacing parameter, r s for different values of the plasma atomic-number. The figure consists of phase-separation 
curves for different plasma composition and each curve has two branches. It is observed that, for the dashed branch 
the parameter, T becomes negative and the potential imaginary. The solid branch, on the other hand, gives rise 
to the real electrostatic potential. The part of the phase-separation curve above the horizontal line indicates the 
presence of the Shukla-Eliasson attractive force and part below the horizontal line corresponds to the ordinary Debye 
screening-potential. It revealed that as the plasma density increases (r s increases) at some point, depending on the 
plasma atomic-number, the attractive potential appears and by further increase in the plasma density (increase in 
r s ) at some point, depending on the plasma atomic-number, the attractive potential changes to the repulsive one. 
Note that the at the critical point, a = oo, where attractive force is about to appear, the force between ions is the 
ordinary Coulomb-force as it is expected. Also at the point a = 1/4, where the attractive force is about to disappear, 



the potential is finite and given in Ref. ( 41|). It is remarked from Fig. 1 that, the starting and ending points for the 



Z = 0,2, 6, 12,26 




FIG. 1. Fig. 1 shows the phase-separation diagram of the Fermi-Dirac plasma without (Z — 0) and with Coulomb correc- 
tion (Z ^ 0) in Log-Log scale. The horizontal line in figure indicates the limit above(below) which the potential is attrac- 
tive (repulsive). The solid/dashed branches show the imaginary /real parts of the phase diagram. The thickness of curves was 
used as a scale to show the variation in the plasma atomic-number. 



„ = (4.3, 4.8, 5.3, 6.3) x 10 25 ™ ' Z = 26 



n„ = 1.3 X 10 25 cm -3 , Z = 0, 1, 2, 3 



rB<t>(r) 

x 10 2 




rB<t>(r) 

Q 




10 2 o.o 



(a) 



(/.) 



FIG. 2. Fig. 2 shows the variation of Shukla-Eliasson attractive potential depth and scaled position with respect to change 
in the plasma number density (plot 2(a)) and atomic number (plot 2(b)) in the low density plasma limit. The thin curves in 
each plot corresponds to the Z — case and the thickness of curves was used as a scale to show the variation in the plasma 
atomic-number and the plasma number-density. 



attractive potential is strongly dependent to the plasma composition and moves to larger densities for larger atomic- 
number. This means that, the attractive force in quantum plasmas starts and ends at higher plasma densities for 
higher atomic-number plasmas. It is also revealed that, density range for which the attractive Shukla-Eliasson force 
exists narrows as the plasma atomic-number increases. It is critical to note that the Shukla-Eliasson force for iron 
composition starts at approximate density of p ~ 75.6675gr/cm 3 which is much higher than the zero-pressure solid 
iron with density around 8gr/cm 3 . This is an indication of the fact that the binding mechanism of Shukla-Eliasson 
force which is due to the electron tunneling and exchange interaction effects is quite different from that of the ordinary 
solids. Figure 2 shows the characteristics of the attractive potential for various plasma parameters. It is remarked 
from Fig. 2(a) that the increase in the plasma number-density leads to the weakening of the attractive force for a 
given plasma number-density (Z — 26 for iron composition). As it was noted in above paragraph, the attractive force 
for iron-composed plasma happens at much larger densities than that of the ordinary solid iron (e. g. see Fig. 1(a)). 
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FIG. 3. Fig. 3 shows the regions and variation of the attractive potential parameters, i.e., the position of the potential 
minimum, f min (Fig. 3(a)), and the binding energy, 4>{f m in) Fig. 3(b), for different values of the atomic-number, Z, in terms 
of the reduced Wigner-Seitz redius, r a . 



It is also remarked from Fig. 1(a) that the ion separation in Shukla-Eliasson quasi-lattice is almost independent of 
the plasma number-density variation. Figure 1(b) shows the variation in potential profile parameters (the binding 
energy and the quasi-lattice parameter) due to change in the plasma composition parameter, Z . It is clearly observed 
that both the binding energy and the lattice parameters are altered significantly due to the change in the plasma 
atomic-number. From Fig. 3(a) it is observed that the equilibrium distance, f min , and the plasma binding energy 
(the depth of the potential minimum, 4>(r m in)) increases strongly as the plasma-density (or the reduced Wigner-Seitz 
radius) increases for all values of the atomic-number, Z . Figure 3(b) further reveals that the plasma binding-energy 
(the depth of the potential minimum, (f>(f min )) decreases strongly as the plasma-density (or the reduced Wigner-Seitz 
radius) increases for all values of the atomic-number, Z. However, it is remarked that the attractive potential occurs 
at higher densities for higher atomic-number values. This can be an indication of the likelihood of presence of heavier 
compounds at the cores of large planets like Jupiter in our solar system. Figure 3 also confirms that there are two 
(higher and lower) cutoff plasma density values for every given atomic-number value which limit the occurrence of the 
Shukla-Eliasson attractive force. It also observed that the attractive force starts closer to the test charge and with 
higher binding-energy values for higher atomic-number of plasmas. The increase of the plasma number-density is also 
observed to decrease the ion separation leading to formation of more compact with higher melting point quasi-lattices. 
Current research highlights features in super-dense degenerate plasma which can be of vital importance in the study 
of planetary sciences. 



CONCLUDING REMARKS 



We used the linear dielectric response of a dense degenerate plasma of electron number-density in the range of 
that in the cores of big planets to show that the Shukla-Eliasson attractive force exists in such density regimes. It 
was found that the nature of newly discovered attractive force in degenerate plasmas is quite different from those 
in ordinary solids and the binding length and strength in this density regime significantly depends on the plasma 
density and composition. It was observed that heavier elements in the plasma give rise to stronger ion-binding. It 
was also remarked that the heavier the plasma ions the denser should the plasma be in order for the plasma binding 
to occur. This research highlights features which can be of fundamental importance in the study of the planetary core 
structures as well as inertial confined laboratory plasmas. 

The author expresses his deepest gratitude to Prof. Padma Kant Shukla for his generosity in exchanging his 
precious ideas. Insightful comments of anonymous referees which have led the paper to its current state are also truly 
appreciated. 
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